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Abstract 



We consider the problem of computing form factors of the massless XXZ Heisen- 
berg spin-1/2 chain in a magnetic field in the (thermodynamic) limit where the 
size M of the chain becomes large. For that purpose, we take the particular 
example of the matrix element of the operator a z between the ground state and 
an excited state with one particle and one hole located at the opposite ends of 
the Fermi interval (umklapp-type term). We exhibit its power-law decrease in 
terms of the size of the chain M, and compute the corresponding exponent and 
amplitude. As a consequence, we show that this form factor is directly related 
to the amplitude of the leading oscillating term in the long-distance asymptotic 
expansion of the correlation function (ofer^ . j.). 
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1 Introduction 

The main purpose of this article is to open a way for the study of the form factors (or in general 
of the matrix elements of local operators) of the massless XXZ Heisenberg chain in a magnetic 
field in the limit where the size of the chain becomes large. For that purpose we will consider 
a particular matrix element of the local spin operator a z between the ground state and an 
excited state with one particle and one hole located at the opposite ends of the Fermi interval 
(umklapp-type term). Although the present article is devoted to this special case, it will be 
clear that the method we propose can be applied to more general cases as well. 

As already mentioned in our recent work [JJ , this particular form factor is of direct interest 
for the computation of the asymptotic behavior of the a z two-point correlation function. Let 
us briefly recall how it appears in this context. In we explained how to compute, from first 
principles, the long-distance asymptotic behavior of some two-point correlation functions: the 
longitudinal spin-spin correlation function in the massless phase of the XXZ spin-1 /2 Heisenberg 
chain and also the density-density correlation function of the quantum one-dimensional Bose 
gas. In particular, for the spin chain, we obtained that the correlation function of the third 
components of spin behaves at large distance m as 



In this formula p F stands for the Fermi momentum and Z represents the value of the dressed 
charge on the Fermi boundary, the precise definition of these quantities being postponed to 
Section [2~2"1 We announced in [1J that the coefficient F a in (|l.ip is related to the thermodynamic 
limit of the properly normalized aforementioned special form factor of the a z operator. We 
prove this statement in this article. 

The fact that the coefficient F a should be related to this form factor is not really surprising. 
It is well known that the leading asymptotic behavior of the correlation functions is defined by 
low-energy excitations in the spectrum of the Hamiltonian. From this assumption, the form of 
the asymptotics (jl.ip was predicted by the Luttinger liquid approach [21 [HJ 01 El [6j [8] and 
the conformal field theory [9j [Til ETUI [T2] together with the analysis of finite size corrections 
[13 HH [13 [TBI HZl HSl [fg HOI I2I]. There the osculating term in the asymptotic formula 
corresponds to the process of a particle jumping from one Fermi boundary to the other. One 
can expect therefore that the numerical coefficient in this term should be somehow related to 
the corresponding form factor. 

However, the precise relationship between the coefficient F a and the form factor under 
consideration has been missing up to now. The first reason is that the methods mentioned 
above hardly give any prediction for the exact value of the coefficient F a : an expression for this 
amplitude was already given in [3(J|, 131] . but only for the case of zero magnetic field, and its 
physical meaning was not clear. The second reason is the absence of sufficient information about 
the thermodynamic limit of form factors in massless models. The matter is that, in distinction 
of massive models in infinite volume for which form factors were successfully calculated in series 
of works (see e.g. [22 1 1231 24, 25J), the thermodynamic limit of the form factors in the massless 
case have a non-trivial behavior with respect to the size of the system. This phenomenon 
was observed, apparently for the first time, in [28] for the model of one dimensional bosons. 
Recently similar results were obtained in [29J for the free fermion limit of the XXZ spin chain. 
This non-trivial behavior of the form factors with respect to the size of the system makes very 
difficult their analysis directly in the infinite volume contrary to the case of massive models. 




corrections. 
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One of the main goal of the present article is to develop a method to resolve this problem. 
We start with determinant representations for form factors for the finite Heisenberg chain 
|36| . Then, acting in the spirit of works [2H\ 125] . we proceed to the thermodynamic limit 
of the specific form factor under consideration. This method allows us to obtain its leading 
(power-law) behavior in terms of the size of the chain and to calculate the corresponding finite 
amplitude. Remarkably, the exponent governing the power-law decrease of the norm squared 
of this (umklapp-type) form factor in terms of the size M of the chain is equal to the critical 
exponent 2Z 2 for the corresponding oscillating term in (II. ip . 

The article is organized as follows. In Section[2l we recall some well known results concerning 
the XXZ chain, its algebraic Bethe ansatz solution and its thermodynamic limit. The main 
result of the article is discussed in Section [3) we define more precisely the form factor that 
we investigate and we write explicitly its relation with the amplitude appearing in we 
also recall the exact expression of this amplitude as computed in [I]. The next sections are 
devoted to the proof of this result. In Section 2] we remind determinant representations for this 
form factor in the finite chain that are suited for our analysis. It enables us, in Section \5\ to 
compute its thermodynamic limit. Several technical results concerning this limit are given in 
the appendices. 



2 The XXZ chain: general results 

The Hamiltonian of the integrable spin-1/2 XXZ chain, with anisotropy parameter A and in 
an external longitudinal magnetic field h, is given by 

H = H<® - hS z , (2.1) 

where 

M 

H {0) = Y, + «+i + A(<4<Cfi - 1)} , (2.2) 

m=l 
1 M 

Sz = ^J2<, [H^,S z }=0. (2.3) 

m=l 

In this expression, as well as in crm V ' z stand for the local spin operators (in the spin-1/2 

representation where they are represented by Pauli matrices) acting non-trivially on the m th 
site of the chain. The quantum space of states is 7i = <8>^ =1 H m where TL m ~ C 2 is the local 
quantum space at site m. We assume periodic boundary conditions and, for simplicity, the 
length of the chain M is chosen to be even. Since the simultaneous reversal of all spins is 
equivalent to a change of sign of the magnetic field, it is enough to consider the case h > 0. 

In the thermodynamic limit (M — ► oo) the model exhibits different regimes depending on 
the values of the anisotropy A and magnetic field h [32]. We focus on the massless regime 
(|A| < 1) and assume that the magnetic field is below its critical value h c [39]. All along this 
article we use the parameterization A = cos £, < £ < ir. 



2.1 The finite XXZ chain in the algebraic Bethe Ansatz framework 

The central object of the algebraic Bethe Ansatz method is the quantum monodromy matrix. 
In the case of the XXZ chain, it is a 2 x 2 matrix, 
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with operator- valued entries A, B, C and D acting on the quantum space of states 7i of 
the chain. These operators depend on a complex parameter A and satisfy a set of quadratic 
commutation relations driven by the six-vertex i?-matrix. 

The Hamiltonian H^> of the XXZ chain is given by a trace identity [33] involving the 
transfer matrix 

T(A) = trT(A) = A{\) + D(\). (2.5) 

The transfer matrices commute for arbitrary values of the spectral parameter, [T(\),T(p,)\ = 0, 
and their common eigenstates coincide with those of the Hamiltonian. 

The solution of the quantum inverse scattering problem enables us to express local spin 
operators in terms of the entries of the monodromy matrix [36} 137] : 



a 



a rfk— 1 



T k -\-iC/2) ■ tr (T(-</2) a a ) ■ T-\-iQ/2). (2.6) 



In the l.h.s. cr^ denote a local spin operator at site k, whereas a a appearing in the r.h.s. should 
be understood as a 2 x 2 Pauli matrix multiplying the 2x2 monodromy matrix. 

It turns out to be convenient to introduce a slightly more general object [26] , the twisted 
transfer matrix 

T K {\) = A(X) + kD(X), (2.7) 

depending on an additional complex parameter k. For a fixed value of k, these twisted transfer 
matrices also commute with each others, [T K (X),T K (/j,)] = 0, and hence possess a common 
eigenbasis. The latter goes to the one of the Hamiltonian (|2.2h in the k — > 1 limit. 

In the framework of the algebraic Bethe Ansatz, an arbitrary quantum state can be obtained 
from the states generated by a multiple action of -B(A) operators on the reference state |0) 
with all spins up (respectively by a multiple action of C(A) operators on the dual reference 
state (0 |). Consider the subspace HW 2 ^ of the space of states 7i with a fixed number N 
of spins down. In this subspace, the eigenstates | ip K ({fi}) ) (respectively the dual eigenstates 
{tpudfJ'}) |) of the twisted transfer matrix T K {v) can be constructed in the form 

N N 

I Mil*}) ) = II B{jJ.j)\ ), ( MilA) I = ( | [J C{ H ) , (2.8) 

3=1 3=1 

where the parameters hi, . . . , hn satisfy the system of twisted Bethe equations 

yMin}) = 0, j = l,...,N. (2.9) 
Here, the function y K is defined as 

N N 

y^lM) = a ( v ) H sinh(/x fc -v-iQ + n d(y) sinh(^ fc - v + (2.10) 

k=l k=l 

a(z^), d(v) being the eigenvalues of the operators A(v) and D(v) on the reference state | 0): 

a{v) =sinh M (^-^), d{v) =sinh M (zy + f). (2.11) 
The corresponding eigenvalue of T K (y) on | ^> K ({//}) ) (or on a dual eigenstate) is 

rMifi) = fr s[nh{ ^- u - i c) + Kd (u) n ^fr-^-*) , (2 . 12) 
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Note that all this construction is also valid at k = 1, in which case we agree upon omitting the 
subscript n in the corresponding quantities. 

Not all the solutions to the system (|2.9p yield eigenstates of the operator T K {v). The detailed 
classification of the pertinent solutions is given in [35J (see also [27] )■ For those solutions 
corresponding to the eigenstates (the so-called admissible solutions), the system of the twisted 
Bethe equations can be recast as 

ap±» ^-H-K) Ki j = l,...,N. (2.13) 
2.2 Thermodynamic limit 

We outline here very briefly the thermodynamic limit of the model, more detailed analysis can 
be found in @Q1 EU [39] . 

The ground state of the infinite chain with fixed magnetization {a z ) can be constructed 
as the M — > oo limit of the finite chain ground state f|2.8[) at k = 1. Hereby the limit of the 
ratio N/M is equal to an average density D whose value is related to the magnetization by 
(a z ) = l-2D. 

In order to construct the ground state, one first takes the logarithm of the system (|2.13|) at 
K = 1 and N and M fixed, 

N 

M Po ( H ) - #(Pj ~ Mfc) = 27rnj, j = 1, . . . , N. (2.14) 

fc=l 

Here -M/2 < rij < M/2, with nj G Z for N odd and n,- € Z + 1/2 for iV even (recall that M 
is assumed even). The functions po(X) and $(A), 

, /sinh(^ + A)\ n . . , /sinh« + A)\ 

po A = i log — ^ -L , i?(A) = i log . .; , 2.15 

\ v smh(f-A)y V smh «- A )y 

are respectively called bare momentum and bare phase (we choose the principal branch of the 
logarithm) . Among all the eigenstates obtained through (|2.14[) , we define the iV-particle ground 
state to be the state with the lowest energy in the sector with N spins down. It is characterized 
by the (half-)integers rij = Ij = j — (N + l)/2 [ID], and the corresponding roots Xj are real 
numbers belonging to some interval [—q, q] called the Fermi zone. 

At the next step the length of the chain M as well as the number of spins down N are sent 
to infinity at fixed limiting value of N/M. In this limit the Bethe roots Xj corresponding to 
the iV-particle ground state condensate, — Xj = 0(M _1 ) (see e.g. J34J), and the system 
(|2.14p turns into a linear integral equation for the spectral density p(A) □. 

A standard way to derive this integral equation is to introduce a ground state counting 
function £(A). Namely, let A(£) satisfy the functional equation 

»( A ®) - h 5X A $ " Afc ) = 27r " ) ' (2 - 16) 

k=l v y 



1 Although there are evidences that the ground state can be described in that way, strictly speaking, to our 
knowledge, this statement has been proved only in the case A < [401 134] . 
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for a given value of £ and with parameters A& corresponding to the iV-particle ground state. 
Evidently, £(A) is the counting function, since X(jj) = Xj and £(Aj) = j/M due to equation 
(|2.14p . The (finite-size) density of the iV-particle ground state is defined by /5(A) = d£/dA. 
From (I2.16p . it is equal to 



1 - 

2tt P (X) = p' (A) - — £ K (A - X k ), (2.17) 

k=l 

with 

K(\)={>'(\) = -— — . (2.18) 

v ' v ' smh(A + <) smh(A - <) 

In the thermodynamic limit, the iV-particle density /5(A) goes to the spectral density p(X) of 
the infinite chain with support on a finite interval [— q, q] . In this limit, discrete sums of piecewise 
continuous functions can be replaced by integrals via the Euler-Maclaurin summation formula 



N N q 

lim i-V/(A fc )= lim ^YfirH-w)) = I f(X)p(X)dX. (2.19) 

k=l k=l _ q 

Thus, replacing the sum over k by an integral in (|2.17|) . we obtain the integral equation for the 
density of the infinite chain, 

pw + ^j K{ - x -^ p^) d ^ = h p '° {X) - (2 - 20) 

The average density D = lim^M^oo N/M is then given by 

N q 

D= lim -L V" 1 = f p(X) dX, (2.21) 



k=l 
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which defines the integration boundary q in (|2.2U|) . 

Observe that the definition of the discrete density p(X) implies that p(X) can also be defined 
as the limiting value 

p(Xj) = lim %±lhliyM = lim 1 (2 . 22) 

J N,M-*oa Xj + i - Xj N,M-*oo M{Xj + i — Xj) 

The method described above allows one to construct the ground state of the infinite XXZ 
chain at fixed magnetization. For the XXZ chain in an external magnetic field, the magneti- 
zation of the ground state is not an independent variable but depends on the magnetic field 
h. Then the boundary of the Fermi zone q in the equation (|2.20p is defined from the condi- 
tion e(q) = (instead of from (|2.21|) ). where e(X) is the dressed energy satisfying an integral 
equation similar to (12.200 : 

q 

e(A) + — / K(X - /j) e{p) dp = h - 2p' (X) sin ( . (2.23) 
2vr J 
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The condition e(q) = provides the positiveness of the energy of any excited state. We see from 
(I2.23P that the Fermi boundary q depends on the anisotropy parameter A and on the magnetic 
field h, i.e. that q is a function of the parameters of the Hamiltonian (|2.ip . We remind that q 
remains finite for non-zero magnetic fields, while it tends to infinity when h — > 0. Throughout 
the present article, q will be kept finite in all the computations. 

Let us finally recall the definitions of the dressed momentum and charge, which are two 
important functions characterizing the ground state. The dressed momentum is closely related 
to the density as 



p(X) =2ir p{n)&n. (2.24) 



Obviously, the thermodynamic limit £(A) of the counting function £(A) can be expressed in 
terms of the dressed momentum according to £(A) = [p(A) + p{qj\/2ir. The dressed charge 
satisfies the integral equation 

q 

Z{\) + ^ / K(X - (j) Z(jj) dfJL = 1, (2.25) 

and can be interpreted in the XXZ model as the intrinsic magnetic moment of the elementary 
excitations [39]. The quantities Z = Z(q) and p F = p(q) = ttD (Fermi momentum) enter the 
asymptotic formula 

The excitations above the ground state can be constructed by adjoining particles and holes 
to the Dirac sea [3H 133 [33 1331 133 133 HS]- In other words, an excited state above the ground 
state corresponds to the replacement of a finite number of the (half-)integers Ij = j — (N + 1)/2 
parameterizing the ground state by some different ones, not belonging to the original sequence. 
In general, this induces a shift of the Bethe parameters {/i} describing this excited state with 
respect to the ground state Bethe roots {A}, and possible appearance of complex roots. 



3 Statement of the result 



We are now in position to be more precise about the result we prove in this article. 

We study a particular form factor of the a z operator, involving a very specific excited state, 
namely the one with one particle and one hole located on the opposite ends of the Fermi zone. 
For definiteness, we consider the case where the spectral parameters of the particle and hole 
are respectively equal to ±g. This means that the (half-)integers which characterize this state 
through (|2.14l) . that we will denote (j = 1, ...,N), coincide with the ground state ones 
(I'. = 1^ for j = 2, . . . , N, but differ at j = 1: I[ = N + 1 - (N + l)/2. There exists another 
way of describing the same eigenstate: by shifting the ground state (half-)integers by 1, i.e. 
by setting I'- = j + 1 — (N + l)/2. In the framework of this approach, we can consider the 
aforementioned excited state as a limit of a twisted ground state. Indeed, similarly to (I2.14|) 
different eigenstates of the twisted transfer matrix are characterized by sets of (half-)integers 
hj in the logarithmic form of the twisted Bethe equations (|2.13p , 

N 

Mpoifij) - ^2^(pj - Mfc) = 27r "-j - *A j = h---,N, (3.1) 
k=i 
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where k = e@ . In the case hj = Ij, we call the corresponding eigenstate the twisted ground state. 
In that case and at f3 = 2ni, the equations (|3.1f) reduce to (12.14|) with rij = Ij = (N+l)/2, 
and, thus, the twisted ground state becomes the excited state described above. 

Our result is the following : 

In the thermodynamic limit, the normalized form factor of the a z operator between the 
ground state | ^({A}) } and the excited state ( ip({fi}) \ defined above is related to the amplitude 
\F a \ 2 (which is a finite number) appearing in the asymptotic formula (jl.ip for the longitudinal 
spin-spin correlation function as 



2 



In other words, it means that the modulus squared of the corresponding form factor behaves as 

\F a \ 2 (J^) for M large. It is remarkable that the exponent 2Z 2 governing the modulus 

squared of the form factor power-law decrease in terms of the size M of the chain is exactly 
equal to the exponent for the power-law behavior of the corresponding oscillating term in the 
correlation function (jl.ip . but there in terms of the distance m separating the two spin opera- 
tors. 

This result will be proven throughout Sections d] and There, the above special form factor 
will be computed and its thermodynamic limit will be taken. 

In order to complete this statement, let us recall the explicit value of the amplitude |-F CT | 2 
obtained in pQ: 

\F,\ 2 = ( 2G(2 '^ SinPf ) 2 - [2vrp(g) sinh(2 9 )] ~** ■ e*"<* • A(P)\, =2m - (3.3) 

Here the notation G(a, x) represents the product of Barnes functions G(a, x) = G(a+x)G(a—x) 
(cf IA.2p . The constants Co and C\ are defined in terms of the dressed charge Z: 

2 J tanh(A-/i) r J tanh(g - A) v 1 

-q -q 



Finally, the coefficient A(f3) is given in terms of a ratio of Fredholm determinants: 



B rs , -, det \l+ o^?C/_„ (w.w 1 )] 



det [/ + &K] 



(3.6) 



In this expression, the integral operator I + K/2n occurring in the denominator acts on the 
interval [— </] with kernel (I2,18p . whereas the integral operator appearing in the numerator 
acts on a contour surrounding [—<?,(/] with kernel 

U_ q {W,W ) - e (3z(w+i() _ e f3+/3z(w-iO ' V-') 
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where 

K K (X) = coth(A + «C)-«coth(A-<), k = e p . (3.8) 

The function z(w) appearing (13. 6|) . (13. 7p . is the m-periodic Cauchy transform of the dressed 
charge Z, 

z(w) = — I coth(A - w) Z(X) dX. (3.9) 

Observe that it has a cut on the interval [— q, q] and satisfies the following properties: 

z(w + i() -z(w-i() = 1-Z(w), (3.10) 
z + (w) — z^(w) = Z(w) , w£[—q,q\. (3.11) 

The equation (|3.10p follows from (12.250 . In (|3.1ip . z± stand for the limiting values of z when 
R is approached from the upper and lower half planes. 

Remark 3.1. Equation (|3.3|) contains the limiting value of A{(3) at (3 = 2iri. We did not set 
(3 = 2-Ki directly in (|3.6p - (|3.8p . as it may cause problems at some particular values of A. For 
example, at A = (C = §, Z(X) = 1), the kernel (13. 7p becomes ill-defined when (3 = 2iri. It 
should therefore be understood in the sense of the limit (3 — > 2in, the latter being well defined. 



4 Special form factor for the finite chain 

The solution of the quantum inverse scattering problem (|2.6H and determinant representations 
for scalar products \38\ [36] give a possibility to express, for the finite chain, arbitrary form 
factors of local spin operators as finite-size determinants |36j. We explain here how to obtain 
such a representation for the special form factor defined above. Then, following the ideas of 
[251 IT], we transform the obtained determinant to a new one, more convenient for the calculation 
of the thermodynamic limit. 

From now on, we will consider eigenstates of the twisted transfer matrix only for k belonging 
to the unit circle, i.e. k = where (3 is pure imaginary. This restriction is not crucial, but 
convenient. The matter is that, due to the involution A^(X) = D(X) (f means Hermitian 
conjugation), the operator e~^/ 2 T K {X) becomes self-adjoint for |k| = 1 and A € R. In such a 
case, the roots of the twisted Bethe equations (|3.ip are real or contain complex conjugated pairs. 
In its turn, the involution B\X) = —(7(A) guarantees the duality of eigenstates (^({a 4 })! = 

(-i)^ K (M))t. 

4.1 Special form factor and scalar product 

Consider the following operator Q m , giving the number of spins down in the first m sites of 
the chain: 

1 m 

Sm=2X>-<). (4.1) 

n=l 

It is easy to see that 

e , Qm = "Q + \_H . ^ , K = e ^ (4.2) 
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and we have 
d 



on ^_ P /32« 
op 



= l-a z rn+ll neZ, (4.3) 

/3=27rin 



where the symbol D m stands for the lattice derivative: D m f(m) = f(m + 1) — f(m). Also, 
note that is a polynomial in k which becomes the identity operator whenever (3 = 2-7rin, 

nez. 

The operator (|4.2j) admits, through the solution of the inverse problem ()2.6p . a simple 
representation in terms of the twisted and standard transfer matrices: 

e PQ m =T™{- |) T~ m ( - |) . (4.4) 

We use this representation and consider the following matrix element of e^® m : 

mM)\e pQm m{\})) = (rp.mwri - 1) T ~ m ( - («) 

where IVKIA})) is the ground state in the iV-particle sector, and (^({^Dl is an eigenstate of 
the twisted transfer matrix T K . Then 

JQmiMiww - ( T "(-rlO*»V 



(MM)\^ m mm = ( "t ( , ) <^(M)IV-({A})>. (4.6) 



Setting v = — ^ in (|2.12j) . we obtain 

(^({/iDle^l^A})) = e im ^=i^^)- p °( A ^](V K ({/i})|V({A})). (4.7) 



Let us now differentiate (|4.7p with respect to /3 at /3 = 2ni. Hereby, at f3 = 2ni, the state 
{tp K ({^})\ becomes an eigenstate of the standard transfer matrix. As it does not coincide with 
the ground state, it is orthogonal to the last one, and thus, 

WM)|Sm|V({A})> = ( e -£f =1 bofe)-Po(A,)] _ i) ^<^({/4)|V>({A}))| • (4-8) 



Using also (|4.3|) we find that the special form factor we want to compute can be obtained in 
terms of the normalized scalar product between the ground state | ?/>({A}) ) and the twisted 
ground state {ip K ({[i}) | given by (|3,ip with hj = If 

l^(M)fe + il^({A}))l 2 fi . 2 m,\ <9 2 cW 



ll^(M)ll 2 -ll^({A})|| 2 

with 



!sin2 (^)-^ )(M ' {A}) 



(4.9) 

/3=2m 



Sn (M ' {A}) " lll^(M)ll-W{A})||J ' (410) 



and where "Pe^ = ^2jLi\po(^j) — Po(Aj)] is the total momentum of the excitation. Deriving 



79]) we have used that (V'ft({^})|V ; ({A})) is real at \k\ = 1. 
In (|4.10p , the parameters fij satisfy the system of twisted Bethe equations (13. ip with -hj = Ij 
and, hence, are functions of (3. In particular, they are taken to be 

H(P)\p = o = X J> j = l,...,N. (4.11) 

Then, in the l.h.s. of (|4.9[) . the set {fij((3 = 2tti)} describes the excited state considered in 
Section [3] with a particle at q and a hole at —q. 
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4.2 Determinant representation 

Thus, in order to study our special form factor, we can use the determinant representations 
for the scalar product between an eigenstate of the twisted transfer matrix with any arbitrary 
state of the form (|2.8j) . Let us recall these representations. 



Proposition 4.1. 138[ \36ij Let m, . . . , /i/y satisfy the system (|2.9p and Ai, . . . , A/v be generic 
complex numbers. Then 

N N 

(oin^)i^(M)) = (^M)in^i)io> 

j'=i i=i 

na=1 ^ a) detArJUMJAUM), (4.12) 



N 

n sinh(/i a - fi b ) sinh(A b - A a ) 

a>b 



where the N x N matrix {A}|{/u}) is defined as 

N 

(tt K )j k ({//}, {A}|{/u}) = a(Afc) Ajfc) sinh(/i a - A fc - <) 

a=l 

AT 

- /cd(Ajfe)t(Afe,/ij) sinh(/i a - A fc + iC), (4.13) 

a=l 

A) = -T7 ~'. S1 ^ C : -r. (4.14) 

^ ' sinh(// - A) smh(/j - A - O v ; 

Note that, due to the Bethe equations (12. 9p . the entries of the matrix Q K ({ft} , {A} | {fi}) 
are not singular at A& = fij. In particular, if Xj = fij for all j = 1, . . . , iV, then we obtain the 
square of the norm of the twisted eigenstate {4>k{{^})\4>k{{iA)) (recall that \k\ = 1), which 
reduces to 

N 

N J] sinh( / u a - fib- i() 



a,b=l 



TV 

j=1 11 sinh(/i a - fi b ) 



a, 6=1 



xdetjvO^, (4.15) 



with 



a=l 



In the case of interest (scalar product of the untwisted ground state | ^({A}) ) with the re- 
twisted one ( ip K ({fi}) |), one has two different representations for the scalar product. On the one 
hand, it is possible to apply equation (|4.12[> when | ^({A}) } in understood as an arbitrary state. 
On the other hand, one can consider | ) as an eigenstate (for k = 1), interpret ( ip K ({fi}) | 
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as an arbitrary state, and apply formula (|4.12|) with A <-> fx. Thus, one can obtain two different, 
but equivalent representations for the scalar product (ip K ({(i}) [^({A})}. We do not give here 
the explicit formulae (which can be easily derived from (|4.12p - (|4.14p ) as these representations 
are not convenient for the calculation of thermodynamic limits. In order to obtain new formulae 
that are appropriate for such goal, one can extract the products of sinh~ 1 (/i^ — A,) from the 
determinant of the matrix f2 K (see [Tl 1281 125] ) . As a result, one obtains two new representations 
containing Predholm determinants of integral operators of the form / + t^Uq^ (w, w'), with 
kernels 



U { e X \w,w') 



-i^r sinh(w - Ha) K K (w - w') - K K (9 - w') 
A = l sinh(u; - A a ) ' V+\w) - kVI\w) ' 



N 



smh(w' - A a ) K K {w - w') - K K (w 



^■j^ sinh(w' — (L a ) V-(w') — kV+(w') 
Here 6 is an arbitrary complex number, K K is given by (|3.8|) . and 



(4.17) 
(4.18) 



v±H = v ± 



w 



N 

n 



sinh(u> — A a ± iQ 
L sinh(w - fi a ± <) ' 



(4.19) 



These integral operators act on a contour surrounding the points {A} (resp.{/x}) but no any 
other singularity of Un^' (w,w'). More precisely, one has the following proposition (see Ap- 
pendix A of [1] for the proof): 



Proposition 4.2. fl^ Let Ai,...,Ajv satisfy the system of untwisted Bethe equations and 
fii, . . . , (an satisfy the system of n-twisted Bethe equations (|2.9j) . Then 



w{a»i^(m))= fi si t^ Aa -j° -n^^)^) 



a, 6=1 



3=1 



1 — K 



vi 1 ie) - kvz\o) 



■ det 



1 



and 



(4.20) 



w«A()fc(w»= n ^t^'ni*'*' 

a, 6=1 



3=1 
1 — K 



V-{0)-kV+{0) 



■ det 



(4.21) 



We would like to emphasize that equations (|4.2(jp and (|4.2ip describe scalar products for 
the finite XXZ chain even though they contain Predholm determinants of integral operators. 
In fact, it is not difficult to see that the integral operators Uq^'^ are of finite rank and hence 
the Fredholm determinants can be reduced to usual determinants of N x N matrices [T]. 

Remark 4.1. We gave here two different representations for the same quantity, because we will 
use both of them. Namely, we will substitute (I4.20p for one scalar product (V ; k({/ x })IV'({A})} 
in (|4.1(jp and (|4.21|) for the other one, despite the evident fact that these scalar products are 
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the same. This way seems to be slightly strange, but it leads directly to equation (|3.3|) , Due 
to the equivalence of the representations (I4.20p and (14.211) one can use, of course, only one of 
them, coming eventually to the same result. The last way, however, requires some additional 
identities, since the equivalence of the representations (I4,20|) and (|4.2ip is non-trivial. 

It remains to apply all these results to our special form factor, which was expressed in 
terms of the normalized scalar product sff through (|4.9[) . Substituting (|4.20p . (|4.21|) . (|4.15l) 



in (|4.10p , and using the Bethe equations (I2.13|) , we obtain a representation for sff as a product 
of three factors: 



S£ } ({/*}, {A}) 



exp Cq I , 



(4.22) 



with 



D 



N 



det 



sinh( / u i - A fc ) 



N 

n 



«V+(A,) 



4vr 2 M 2 p(A i )p K (/x. 



•3) 



(4.23) 



and 



C, 



0,N 



N 



a,6=l 



sinh(A a - /ifc - sinh(^ b - A a - Q 
sinh(A a - X b - i() sinh(^ a - fj, b - i() ' 



(4.24) 



A 



N 



(k - I) 2 det 



I + ^\w,w') det J+^C/^K^) 



27ri^6» 2 



(y-i(^) _ wKT 1 ^!)) (V-(9 2 ) - kV + (9 2 )) ■ det N • det^ 8$ 



(4.25) 



Recall that the representations (|4.20p . (|4.2ip contain an arbitrary complex parameter 9. Here 
we have used two different parameters 9\ and 9 2 for each of the scalar products. In the following, 
we shall set 9\ = —q and 8 2 = q- 



5 Thermodynamic limit of the special form factor 

We compute now the large size behavior (N, M — > oo, N/M — > D) of the special form factor 
represented in (|4.9p , (|4.22p . We need first to characterize how the parameters \ij of the k- twisted 
ground state Bethe Ansatz equations behave in this limit. More precisely, how they are shifted 
with respect to the parameters Aj of the untwisted ground state. This will enable us to take the 
thermodynamic limit in the expression (|4.22p for the normalized scalar product. Whereas it is 
quite straightforward to evaluate the limit of the factors A$ and (see Subsection 15. 2p , 

we will see in Subsection [O] that the factor D$ (|4T23j) requires more attention. 

5.1 Twisted ground state and shift functions 

In order to study the K-twisted ground state |V ; k({/^})) in the thermodynamic limit, it is 
convenient to introduce, similarly as for the untwisted case (see (I2.16P ). a twisted counting 
function £ K (^) satisfying £ K (/Uj) = j/M. Its explicit expression is given by 

k=l 
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with pq and # given by (|2. 15 j) . The corresponding discrete iV-particle density p K {p) = ^ 
coincides with the function defined in (|4. 161) . It also coincides, at (3 = 0, with the iV-particle 
density of the ground state /5(A). It is easy to check that both functions p K {X) (for k / 1) and 
p(A) go to the same function p(X) (|2,22[) in the thermodynamic limit. 

Using the counting functions for the ground state (|2.16p and twisted ground state (|5.ip . 
one can define shift functions describing the displacement of the spectral parameters {p} with 
respect to ground state ones {A}. One can actually introduce two shift functions: 

F(Xj) = M - aXj)] , F K (pj) = M [Uh) ~ U*s)] • (5-2) 

Both have the same thermodynamic limit which we, from now on, denote by F(X). One can 
easily verify that 

F(X j )= lim ^~ X { , (5.3) 

J N,M-^°° Xj + i - Xj 

The subtraction of the equation (|2.14p with rij = Ij (Bethe equation for the ground state) from 
(|3.1|) (Bethe equation for the twisted ground state, corresponding to the same number hj = Ij) 
and a replacement of finite differences by derivatives and of sums by integrals as in (|2.19p , leads 
to the integral equation for the shift function F: 

q 

F{x) + hf K{x ~ M)F(/x) d/i = 4h- (5 - 4) 

-q 

Comparing this equation with the one for the dressed charge (|2.25j) . one obtains 

Thence, for M large enough, the spectral parameters pj are shifted with respect to the ground 
state ones Xj by a quantity ej = pj — Xj which can be computed for large M 



27riMp(Xj) 



+ 0(M- 2 ), j = l,...,N. (5.6) 



5.2 Thermodynamic limit of and Cq K 1, 



(k) t (k) 

We first compute the limits of Cq L and A N . These limits are quite simple to calculate since 



it is enough to take into account only the leading order of the shift (|5.6|) . 

We have already seen that, in the thermodynamic limit N, M — > oo, discreet sums over the 
parameters Xj can be replaced by integrals via (|2.19j) . It is also easy to see that 



N / 1 \ 
lim TT 1 + —f(Xj) = exp 

3=1 V 7 



lim det at 

N,M-*oo 



6 jk + jjV(^j,h) 




(5.7) 



(5.8) 



provided /(A) and V(X, p) are piecewise continuous on [— q, q\. In (|5.8|) . det[I + Vp] refers to 
the Fredholm determinant of the integral operator / + Vp acting on the interval [— q, q] with 
kernel V(X,p)p(p). 
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( K \ 

The thermodynamic limit of Cq ^ is straightforward to obtain using (15, 6\\ : 



lim C {R) = lim V log sinh ( A i ~ X k ~ g fc ~ K) sinh(Aj - X k + - tQ 
n,m^oo °> N n,m^co Z—^ sinh(Aj — Xk — i() sinh(Aj — X k + ej — e k — i() 



P 2 



lim ej e k d\.d\, logsinh(Aj - X k - i() = — T C , (5.9) 



j,k=l 



where Co is given by (|3.4p . 

Let us now consider the different factors in f|4.25[> . As the limit of /5(A) and p K (A) coin- 
cides with the ground state density p(X), both determinants detjvS^ go to the Fredholm 
determinant of the operator / + K/2tt: 

lim det^v 9^' M) = det [I + K/2k] . (5.10) 

N,M^oo ■> 

A direct application of (15. 7|) leads to 



lim F±M = e ~/M™±0. (5.11) 



' — >oo 

Similarly. 



lim TT ; ^ = e P z ( w ^ id uniformly away from [— o, ol. (5.12) 

N,M^oo A - 1 - sinhfu; — u a ) 
a=l 

Substituting f)5.11[) and (|5.12p into (I4.17P we have that, at 9\ = —q, 

lim J_u^(w,w') = ^-U^(w,w'), (5.13) 

where the kernel U_g(w, w') is given by (13. 7p . It is also easy to see that, when 02 = q and f3 is 
purely imaginary, 

lim — U^(w,w') = (— U (X A\-w,-w'), (5.14) 
where f means hermitian conjugation. A simple algebra based on (|3.10p leads to 

lim A P = (^Jrl (5.15) 
iV,M-+oo w ^sinh^ / 

where *4(/3) is given by (13. 6p . 



5.3 Thermodynamic limit of the Cauchy determinant 

As already mentioned, the thermodynamic limit of the factor (|4.23p is the most com- 

plicated to obtain. Although the final result for this quantity can be expressed in terms of 
the limiting value of the shift function (|5.5p . the approximation (|5.6|) is not enough for the 
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(k) 

intermediate analysis. The reason is that D N contains the square of a Cauchy determinant 
which has poles at fj,j = A&. These singularities are compensated by the zeros of 

km- o 

However, this compensation takes place if and only if the parameters jXj exactly solve the set 
of twisted Bethe equations, which means that it is not enough to consider only the first order 
of the shift function as in (|5.6|) . Therefore, to compute the limit of D N , we should work with 
the exact formulae for the iV-particle shift functions F and F K (|5,2p defined in terms of the 
two counting functions and (, k (uj) (|2.16|) . (|5.1|) . their limiting values being only taken at 
the end of the computation. 

In order to re- write (|4.23p in terms of the shift functions F and F K , we also introduce the 
functions 

. sinh(A — a) A ,, , sinhfA — u) ,„ ,_. 

^ (A ' M) = ?m / ' ^ k(a ' m) = rm rr? (5 - 17) 

6(A) -60") 6«(A)-6«(m) 

which fulfill 

^(A.A) = 6'(A) =p(A), ^(A,A) = &(A) = p«(A). (5.18) 
Note that and <p K have the same thermodynamic limit denoted by <p. 
Proposition 5.1. The product D$ f|4.23j) admits the following factorization: 

< = H[F] ■ H[-F K ) ■ m • *[<?*] ■ ft ^4#^TT - ^ 
The functions H[F] and $[<p] are defined as 

+ , 20) 



AT 

EI [<p(Ai,Aj|,)^(Mfc,Mj)] 



= , (5-21) 

n <p(pj,h) 
j,k=i 

and the function H[—F K ] (respectively $>[<p K ]) can be obtained from (15 .201) (resp. (I5.2ip ) by the 
replacement F(X) — > —F K (fi) (resp. tp — > <£ K/ ). 



Proof — Using the definition (I2.16|) . (|5.ip of the counting functions, one can easily see that 

kV + (uj) _ e 2xiM[i{u>)-i t {u>)] j ^^2) 
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which leads to 

N 



n 



«V+(A 3 -) 



KV + (flj) 
N 



5=1 
N 



-27riAf[|(Mj)-l«(w)] 



JQ- ^ e a»r*JWrK w O*i)-fi ( .CA J )] _ ^ e -2 J riMK(A» J -)-«Aj)] _ x 

5=1 
N 

\\ [4sm(vrF K (^)) sin( 7 r#(A i )) 



5=1 



(5.23) 



Here we have used that = £«(^j) = j/M, and the identity 



AT 



(5.24) 



5=1 



which is a consequence of (I2.16D . 115. ip . 

The remaining part of the proof is quite trivial. It is based on the expression for the Cauchy 
determinant in terms of double products 



det 



JV 



1 



Uj - v k 



JV 

n - u k ){v k - Vj) 

j>k 

JV 

n (uj-v k ) 
j,k=i 



(5.25) 



where Uj, Vj are arbitrary complex numbers. Substituting the explicit expressions for all the 
factors into (|5.19j) . and using again = C«(a*5') = j/M together with the definition of F 

and F K (j5.2|) . one easily reproduces the original representation (|4.23p for D 



(ic) 



□ 



Since the functions (p(\, /i) and <p K (\,[i) are not singular at A,// € [— q,q], the thermody- 
namic limits of the last three factors in the r.h.s. of (15. 19H can easily be calculated similarly as 
in the previous subsection. Using (|5.6p we find 



JV 



lim <1> •'!» ,:, • 1) 



N,M^oo 



exp 



4vr 2 



J Z(A) Z(m) d\d^ log p(A, m) dA d/i I . (5.26) 



It remains to compute the limiting values of H[F], H[—F K ]. To do this, it is convenient to 
split H[F] into two factors: H[F] = H^F] ■ H 2 [F], with 



Hi[F] = exp 



JV 



52m[H*j)-H*N-j+i)] 

5=1 



JV 



(5.27) 
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H 2 [F] = exp 



N 



J2^j)[F(Xj)-F(X N ^ +1 ) 



j,k = l 



sin[vrF(A j )] 



(5.28) 



V'(z) being the logarithmic derivative of the r-function (c/ Appendix I A. 1|) . 
Proposition 5.2. XYie thermodynamic limit of the factor H\[F] (|5.27p is given by 



lim log ili [Fl 



(5.29) 



Proof — Due to the smoothness of F and £, there is an e, < e < 1, such that, for M large 
enough, 



F(A,)-F(A fe ) 



j -k 



1 

M 



F(A,)-F(A A 



e(Ai)-^(A fe ) 



< e , uniformy in A,- , A^ . 



(5.30) 



The logarithm of the double product in the r.h.s. of (|5,27l) can therefore be expanded into its 
Taylor integral series of order 2 



N 



log n i + 

j>k V 



F(Xj) - F (X k ) 
j - k 



M 3 





2 

= E 

n=l 


N 


1 


E 


/ dt 




/ ( 



N 



- ( _ ir+ i - FjX^-FjX,) 

nM n | (Ai) _| (Afc) 



J>fc 



F(A,-) - F(X h 



t 



m e(Ai)-€(A fc ) 



(5.31) 



The expressions under the double sum over j and /c being not singular at Xj, £ [— g, g], the 
thermodynamic limit of the integral remaining terms in (|5.31|) vanish due to the M -3 pre- 
factor. The contribution issued from the n = 1 term cancels with the pre-factor in (|5.27|) since 



E F(A, 'lf At) = E 7^t = E *) - ■ 

7,fc=l 



(5.32) 



Thus, the limit of log .Hi reduces to the one of the term n = 2 in the expansion (|5.3ip . which 
ends the proof. □ 



Proposition 5.3. The thermodynamic limit of the factor Hz[F] (|5.28|) is given by 

(F 2 l-q) + F 2 (q)\ , . . F(-q)-F(q) ,-. , „ \ 

H 2 [F] ■ N 2 )=G(l + F(-q)) G(l - F(q)) • e 2 (l-l°g^) 



x exp 



F 2 (A)-F 2 (g) F 2 (A)-F 2 (- 



e(A)-e(?) c(A)-e(-?) 



p(A)dAV, (5.33) 



where G(z) is the Barnes function (1A.7|) . 
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Proof — The product over k in (|5.28|) can be expressed in terms of T-functions as 



iiA M ^) iir(7 + ^ , (A i ))r(i-JF i (A JV _ i+1 )) 



(5.34) 



This means that i?2[-^1 can itself be factorized into two factors, 
H 2 [F] = 4 1) [F].H^[F], 

with 



(5.35) 



(5.36) 



H2 [F] being obtained from (15.360 by the replacement F(Xj) — > — F(Ajv-i+i)- 

We first consider H^[F]. Let F_ = F(| -1 (0)). Evidently, F_ -> in the thermody- 

namic limit. We use the fact that the shift function F(Xj) is bounded uniformly for any value 
of N. Thus we can use the following bound |-F(A,)| < n for some n G N and j = 1, . 
Then we represent [F] as 



,2V. 



i=i *C7 + ^-y i=1 
Using ([ATO]) . (|R4|) we find 



r(i + nAi)) 



r(j + ^(A,)) 



AT 

n 



r (j) e ^(j)^- g(2V + 1) G(l + F_) 



f^T{j + F.) G(2V + 1 + F_) 

and thus, applying the asymptotic formulae ()A.3[) . ()A.9|) we arrive at 



(5.38) 



lim iV" TT 

N,M-*oc r(j + F. 



G(l + F(-g))-e 



S^i(l-log27r) 



(5.39) 



The limit of the second finite product in (|5.37)) evidently is equal to 1 due to F(Xj) — F_ = 
0(M _1 ) for j = 1, ... , n — 1. Finally, in order to compute the limit of the last product in 
(|5.37p we can use Taylor integral representation for the logarithms of T-functions (recall that 
\F(X j )\<n<j) 



" T(j + fUf^- P -WW 

; \ = exp < 

fl r(i + F(A,)) 



N 



j=n 



1 



N 



dt(l-t) 2 ^(F{X 3 )t + j)F\X 3 )-4,"{F_t + 3)F* 



(5.40) 



J=n 
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The limit of the first sum in (|5,4U|) is obtained by applying Lemma IB. 21 proved in Appendix 
IBl Indeed, as the counting function is invertible and satisfies = j/M, we write F(Xj) = 

^(£ _1 (if))- Then, using (TB~9l) . we find 



N 



lim ^'(j)[F 2 _-F 2 (X 3 )] 



N,M^oo 



"F^c'^-F'irHo)) 



dx. 



(5.41) 



o 



Recall that D = lirnjvM— *oo N/M is the average density, so that, after the change of variables 
£ (x) = A, we obtain 



N H 

N ^JZ^mF*-F\X ] )]=-j 

i =1 -a 



F\X)-F\-q) 

e(A)-c(-9) 



p{X)dX. 



(5.42) 



Finally, the same transformation on the function F allows us to apply, uniformly in t, 
Lemma lB.31 to the second sum in (I5.40p . hence proving that its thermodynamic limit is zero. 
Substituting these results into (|5.37|) . we obtain that 



lim 

N,M— »co 



H^[F].N^ 



G(l + F(-q)) -exp 



F(-q) 



(l-log2vr)-- 



1 f F 2 (X)- F 2 (-q) 



£(A) -£(-?) 



p(X) dX 



. (5.43) 



(2) A 

The thermodynamic limit of [F] can be deduced from the previous calculation due to 
the fact that F(X N - j+1 ) = ^r^T ~ &))■ D 

The limit of H[—F K ] is obtained in the same spirit as above. Since the limits of F, p and £ 
coincide with those of their K-deformed analogs, the result is given by (|5.29p and f)5.33|) provided 
that one makes the replacement F — > —F. 

This enables us to obtain the thermodynamic limit of D^': 

(re) 

Proposition 5.4. The product D N admits the following thermodynamic limit: 



lim 



fiz 2 



fiz 2 



■Ci 



N,M—>oo 

where C\ is given by fl3.5H . We remind that G(a, z) = G{a — z)G{a + z) and Z = Z(q). 

Proof — Using the expression (|5,5p of F and the symmetry properties of Z, p and £, 

Z(A) = Z(-A), p(A)=p(-A), and £(A) = D - £(-A), 
we obtain from Propositions 15.2 1 15.31 that 



(5.44) 



lim 

N,M—>oc 



G 2 (hH) exp 



£_ f f z{X)-z(p) 



8tt 2 j ve(A)-e(M) 



p(X)p(p) dX dp 



x exp 



B 2 } B 2 } Z 2 (X) - Z 2 

£2 J z(x)z(p) dxd.hg^x, p) dXdp - ^ J -^J—^ p (x)dx 



(5.45) 
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After some simple algebra the combination of integrals in (|5.45|) can be reduced to the constant 
C\ (13. 5h . and the result takes the announced form. □ 



5.4 Results for the normalized scalar product and special form factor 

Putting together the individual limits (|5.44p . (|5.9p and (|5.15p yields 



lim ( S { $ • 



2 sinh 







2 r (n §Z^ 

a ^ V ' 2ni, 



P 2 z 2 _ ^ 



[p(q) sinh 2q\ 2n ' 2 e ^ 



(Ci-Co) 



(5.46) 



The limits being uniform in (3 and its derivatives, we take the second (3 derivative at (5 = 2ni 
and obtain 



lim 

N,M— »oo 



2 rt(re) 



diS 



/3 a N 



/3=2ni 



■ M 



2Z2 \ 1 (G%Z) 



2 e 



C\— Co 



y2Z 2 



/3=2ni 



(5.47) 



2\ ttZ J S i n h2g] 

We still need to compute the excitation momentum V ex in (|4.9|) to finish our proof. We 
have 

(5.48) 

AT," ' " x '' ' " ' • 1 ' ' " 



N 1 

lim Vex = lim V[po(Atj) -Po(Aj)] = / Po(^)^(A)dA. 

i =1 -o 



The integral equations (|2.20p . (|2.25p for the density and the dressed charge then lead to 



J p' (\) Z(A) dA = 2vr J p(A) dA = 2p, 



(5.49) 



due to (|2.24p . Substituting (15.491) and (|5.47|> into (f4T9|) and comparing the result with the 
amplitude (|3.3p . we finally arrive at (|3.2|) . 



Conclusion 

In the present article, we have initiated the development of a method to study form factors in 
the massless regime of the XXZ Heisenberg spin chain. In particular we have given explicitly the 
thermodynamic limit of a special form factor in the XXZ chain. This form factor corresponds to 
the matrix element of the operator a z between the ground state and an excited state containing 
one particle and one hole on the different ends of the Fermi zone. It is clear however that our 
method can be applied to other form factors as well, where excited states contain arbitrary 
number of particles and holes in arbitrary positions. This method can also be applied to other 
integrable models, for example to the system of one-dimensional bosons. 

One of the unsolved problems is to prove that at zero magnetic field our result coincides 
with the amplitude predicted in [301 131] . The limit h — > corresponds to the limit q — > oo 
in (|5,46p . It is easy to see that the constants Co, C\, A, p(q) sinh 2q, being finite for finite 
q, become divergent at q — > oo. We were able to prove that the total combination giving F a 
remains finite for h = 0, but we did not obtain a simple expression for its value in this limit. 
We succeeded nevertheless to compute explicitly this quantity in the vicinity of the free fermion 
point (C = f ) up to the second order in e = ^ — |. Our computation confirms the result [301131"] 
up to that order. 
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A ^-function and Barnes function 

We recall in this appendix the definitions of the ^-function and the Barnes function, as well as 
several standard formulae that are useful for our study. 

A.l The V>function and its derivatives 

The ^-function is defined as the logarithmic derivative ip(z) = — logT(z) of the T-function. 

az 

The multiplication property of the r-function implies that, 

^ n \z + l)-^\z)= y —J—. (A.l) 



z 



It follows from (|A.1D that 

1 _ (- 
(k + a) n+1 n\ 



E7i^TT = ^ L (* , " ) ( iV + °»-''''" , (°))- < A ' 2 > 



fc=0 

When z — > oo with — ir < arg(z) < ir, one has 

, ( ^ logi -l + (JL). ,c, ( ^ (-y-H-^ +Q (JL). (A . 3) 

For n > 1, the n th -derivative of the ^-function admits the following integral representation: 

^(z) = -i xZ ~ ll ° gnX d X . (A.4) 
J 1-x 



The latter implies in particular that, for z > 0, 

(-l) n "V (n) (2) > 0, n>l. (A.5) 

It is also easy to see that 
l 

i//(z) - ~ = J x z ~ x -l) dx>0, z>0. (A.6) 
o 
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A. 2 The Barnes function 

The Barnes function G(z) is the unique solution of 

d 3 

G(l + z) =T(z)G(z), with G(l) = 1 and —,logG(z)>0, z > 0. (A.7) 

az 6 

It has the following integral representation 

2 

log G(l + z) = Z<yl ~ ^ + | log 2-K + j x ifj(x) da, > -1. (A.8) 

o 

The Barnes function has the following asymptotic behavior when z — ► oo, —tt < arg(z) < it: 

log G(l + *) = - 1) log z - ^- + | log 2vr + C'(-l) + O (i) , (A.9) 
Due to fA"77l) . one has 



JV-l 



log r(A; + a) = log G(A r + a) - log G(a). (A.10) 

fc=0 

B Sums with logarithmic derivatives of T-function 

In this appendix, we explain how to compute the thermodynamic limit of some sums involving 
the -^-function and its derivatives. 

B.l Finite sums involving ^-function 

Lemma B.l. For n > 0, a / —1, —2, . . . , and an arbitrary complex a, one has the following 
identity 



N 1 

k=l 



N - 1 p a(k+l) 

^){N+a)e^ N ^-^ n \a)e a -(-lYn\ £ + 



fc=l 

Proof — Denote the l.h.s. of (jB.ip as /(a), 



(B.l) 



AT 

/(a) =J]V (n) (fc + a)e afc . (B.2) 
fc=l 

Shifting k by A; + 1 and using (|A.ip . we obtain 

JV-l 

f(a)=e a Y4 n \k + a + l)e ak 

k=0 

N—l JV-l a (fe+l) 

fc=0 fc=o v ' 

AT— 1 a(fc+l) 

= e a /(a) + e a V (n) (a) - e a[ - N+1 ^ {n \N + a) + (-l) n n! ^ ™. (B.3) 

fc=o ( /c + a )" 

Then (jB.ip follows immediately. □ 
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Taking derivatives of (|B. lj) with respect to a at a = one obtains formulae for the sums 
of the type ^ n \k + a)k p . Let us give explicitly several of them. 
One has for p = 0, n = 0, 



N 



^ 4>(k + a) = (N + a)tp(N + a) - a^{a) - N, (B.4) 



k=l 

and for p = 0, n > 1, 

v 



^V (n) (^ + a) = (iV + a)V' (n) (A r + a)-aV' {n) (a + l) + n(V' (n ~ 1) (^ + a)-V' (n " 1) (a + l)). (B.5) 



fc=i 

One also has for p = 1, n = 1, 

■(^(iV + a)-V(a + l)) +^^, (B.6) 

and for p = 1, n > 2, 

v 



E V , ,//, \ (N + a)(N + 1 — a) . a(a — 1) ,/_\, ,. 

fe^(fc + a) = ^ ^- ^//(iV + a) + -i— — -ip ( n) (a + 1) 
At Zi 

k=l 

2a-l/,„ \ N-l 



£ W (n) (k + a)= (N + aKN + l-a) ^ {N + a) + «(a -l) ^(n) (fl + 1} 

k=l 

_ n(2a-l) ^ (n _ x) (jV+a) _ ^ + ^ _ "(2_ 1 l)^(„-2) (JV + fl) _^(n-2) (a+1) ^ 



(B.7) 



B.2 More complicated sums 

We now consider sums involving in addition some regular function(s). 
Lemma B.2. Let f E C 1 ([0, a]) for a > 0. Let us consider the sum 

N 



ri/]=E[/(i)-/(°)]^ (B-8) 



j(JV,Af) ^ 

fc=l 

In the limit N, M — > oo ; iV/M — > D with D G [0, a[, it tends to 



S { N ' M) [f]^j m - f{0) dt. (B.9) 
o 

Proof — We have 

N N ^ 

sr%i = e H - 1 k)^Jf'( ! M) dt+ T.jiJ /'(f) *■ 

fc=l V 7 n V 7 fc=l n V 7 
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It is easy to see that the first sum vanishes in the limit considered. Indeed, setting 
su P[o,a] l/l> we S et using (]B.6|) . (|A.6|> and (1A.3|) . 



A' 



fc=l 



kj M J J \M 
o 



< 



Ml 

M 



N 



fc=i 



7 11 + l)V'(iV) + ^(iV) _ ^(i) _ at - 1] = M [log iV + O(l)] . (B.ll) 



2M 



2M 



Thus, the limit of s[ N ' M ^ [/] reduces to one of the second term in (jB.lOp . which is quite straight- 
forward due to Euler-Maclaurin summation formula 



lim s[ N ' M) [f]= lim / f'(^-)dt = I dy I f'(yt)dt 



D 1 



D 



N 



k=l 



/(y)-/(o) 



dy, 



o o 



and Lemma lB.21 is proved. 



(B.12) 



□ 



Lemma B.3. Let F, f € C 1 ([0, a}), with a > 0. Let n £ N be such that ||-F|| = sup^^ \F\ < n. 
Then, for m > 1, the sum 



N r , 

S^ M) [f,F] = £ /(^> (m) (^) + *) -/(0)^ m) (^(0) + ^) 



k=n 



(B.13) 



vanishes in the limit N,M — > oo, N/M —* Z? D £ [0, a]. 



Proof — We have 



k=n 



i 



. (B.14) 



Due to CO]), (HD (-l) m -V (m) (^) is 

positive monotonically decreasing function for z > 0. 

Then one has the following estimate 



\S^ M) [f,F]\< 



F'\\ 11/11 + 11/'! 



1/ 



^ fc [(-l) m ^ (m+1) (fc " ll^ll) + (-l) m_1 ^ {m) (A: - 

(B.15) 



k=n 



The last sum is readily computed by (1B.7|) . It is easy to see that it is a O(logiV) at most. □ 
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